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Neutrinoless double beta decay pops up almost in any extension of the standard model. It is perhaps the 
only process, which can unambiguously determine whether the massive neutrinos are Majorana or Dirac type 
particles. In addition from the lifetime of this decay, combined with sufficient knowledge of the relevant nuclear 
matrix elements, one can set a constraint involving the neutrino masses. Furthemore, if one incorporates the 
recent results of the neutrino oscillation experiments, one can determine or set a stringent limit on the neutrino 
mass scale. In addition one may obtain usefuU information regarding the presence of right handed currents and 
the right handed neutrino mass scale. One can also constrain the parameters of supersymmetry and, in particular, 
set limits in of R-parity violating couplings as well as get information about extra dimensions. 



1. THE INTERMEDIATE NEUTRINO 
MECHANISM 

Neutrinoless double beta decay occurs when- 
ever ordinary beta decay is forbidden due to en- 
ergy conservation or retarded due to angular mo- 
mentum mismatch, but the nucleus, which is two 
units of charge away, is accessible via a order 
week interaction P- [21 . It is a process known for 
almost 70 years, which has been searched for, 
but not seen yet. It is still of great theoreti- 
cal and experimental interest since perhaps it is 
the only process, which can unambiguously deter- 
mine whether massive neutrinos can be Majorana 
or Dirac type particles. It can occur only if the 
mass eigenstates are Majorana particles (the par- 
ticle coincides with its antiparticle). 
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where is a 3 x 3 left-handed (isotriplet) Ma- 
jorana mass matrix, niD is a Dirac 3 x N mass 
matrix and TOat is the right-handed (isosinglet) 
N X N Majorana mass matrix. In models moti- 
vated by GUT's TV = 3, TOj, = 0, mo is of the 
order of the up quark mass matrix and toat is very 
heavy (> 10^° GeV) so that hght neutrinos with 
mass of order TOD(TO„)~^m^ can occur. On the 
other hand in R-parity violating supersymmetry 
only nil, is non zero. Anyway after diagonaliz- 



ing the above matrix the weak eigenstates are 
related to the mass eigenstates via the equation: 

In the above notation J/*-^^^ — Umns is the usual 
Maki, Nakagawa, Sakata charged leptonic matrix, 
which for convenience sometimes will merely be 
indicated by U . 

1.1. Left-handed currents only 

In the presence of only left handed currents 
one has a contribution to neutrinoless double beta 
decay arising from the diagram shown in fig. ^ 
One encounters two lepton violating parameters 
■qu and 77^ given by jlj: 
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with nij (Mj) the light (heavy) neutrino masses 



and ^, 



[^1 = e 



the Majorana phases of 



the corresponding mass eigenstates. The separa- 
tion of these particle physics parameters from the 
nuclear physics holds, if the neutrino is very light 
{rrij << TOe) or much heavier than the proton 
{Mj » nip). 
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Figure 1. neutrinoless double beta decay in left 
handed theories for light (a) and heavy (b) inter- 
mediate neutrinos. 
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1.2. The leptonic right-left interference 
contribution 

In the presence of right-handed currents one 
can have a contribution via light intermediate 
neutrinos, in which the corresponding amplitude 
does not explicitly depend on the neutrino mass 
(see Fig. . We encounter two lepton violating 
parameters: 

3 

V = emL,X = nr^RL.VRL = ^ (^^ef ^C^ij'^) 0(4) 

i 

where k = m|/m|j, , e = tanC,. mL,mfi are 
the gauge boson masses and C the gauge boson 
mixing angle. 

1.3. Right-handed neutrino mass term 

In the presence of right-handed currents one 
can have additional mass dependent terms (see 
O. The corresponding lepton violating parameter 
is given by: 

< = + + coe.)f;^ iUl^f ^ (5) 
3 ■' 

1.4. Analysis of the data in terms of light 
neutrinos 

In the presence of right handed currents the 
neutrinoless double beta decay lifetime is given 

n, m by: 
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Figure 2. Diagrams leading to neutrinoless 
double beta decay with amplitude not explicitly 
dependent on the neutrino mass. 
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Figure 3. The heavy neutrino contribution due 
to the right handed current 
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Table 1 

Summary of the present experimental results. 



Isotope 


(y) 


Refs 


«Ca 


> 9.5 X 10^^(76%) 


9 


76Ge 


> 1.9 X 1025(90%) 


m 




> 1.6 X 1025(90%) 


m 


82 gg 


> 2.7 X 1022(68%) 


[12 


lOOMo 


> 5.5 X 1022(90%) 


m 




> 7 X 1022(90%) 
(3.52 ±0.11) X 10-4 
(geochemical) 
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128,130rpg 
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i28Te 


> 7.7 X 1024(90%) 


Hi 


i30Te 


> 1.4 X 1023(90%) 




i36xe 


> 4.4 X 1023(90%) 


ini 




> 1.2 X 1021(90%) < 3(90%) 
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\Xr\^ + |77|2 + Cx^\XHcos{iPi - 1P2)] 
+ \Xr\^ + ReiCoxXXn + Co^r]XR)],i6) 
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where the subscript H indicates heavy particle 
(neutrino), tpi is the relative phase between Xl 
and A and V'2 is the relative phase between Xl 
and 77, while Ci depend on the energy and nu- 
clear physics The ellipses {...} indicate contribu- 
tions arising from other particles such as interme- 
diate SUSY particles, unusual particles predicted 
by superstring models, unusual messengers in ex- 
tra dimensions, exotic Higgs scalars etc. 

As we have already mentioned the experiment 
(see table ^ imposes only one constraint among 
the lepton violating parameters. Normally one 
assumes that only one lepton violating parameter 
dominates and proceeds extracting a limit on that 
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Figure 4. The lightest neutrino mass obtained 
from OvfiP decay and neutrino oscillation data for 
normal (top) and inverted (bottom) hierarchy. 



one (see e.g. Table |31 below). In the most favored 
scenario one assumes that the light neutrino mass 
mechanism dominates. One can then employ the 
the neutrino oscillation data to get the con- 
straints of Fig. ^ Clearly OiyP/3 decay can differ- 
entiate between the hierarchies. In fact the de- 
generate scenario supports the recent claims that 
Oz//3/3 decay may even have been seen ^| . In the 
inverted hierarchy, 0^(3(3 can even set a limit on 
the lightest neutrino mass mi. 

A more reasonable procedure would be to con- 
sider multi-dimensional exclusion plots as, e.g., 
that involving < vii, > and ry, shown in Fig. 
obtained with nuclear matrix elements found in 
the literature [7j. From this plot one extracts the 
standard limits as the intersection of the exclusion 



Figure 5. The excluded region in the space of 
the lepton violating parameters < > and 77 
entering neutrinoless double beta decay. 



Figure 7. The value of e extracted from Oi//3/3 
decay as a function of the average right-handed 
neutrino mass 
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Figure 6. The constraint on the i] A parameters 
extracted from Oi^/3/3 decay. 



Figure 8. The same as in Fig. [71 for the right 
handed boson mass. 



curves from the axes. One may use neutrinoless 
double beta decay to extract information on the 
parameters of the right handed interaction k and 
e from the parameters A and i] (see Fig. (S)). from 
the definitions one finds: 

- = - (7) 
77 rj 

Furthermore Eq. Elcan be cast in the form [P,|ni: 
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(8) 



+ £2 + coeK) < — ^ >= ,7«, Co 

Where <^>=E,^ (t/-)^ S, ^. 

Using Eqs O and © one can extract limits on 
mwR and e. Thus, e.g., using: 

A ^ M X 10-^77 = ^ X 10-«,77« ^ 0.8 X 10-^ 



we obtain the curves shown in Figs and [S] as 
a function of the average heavy neutrino mass. 
These should be compared with the values 

m^R > 715GeV , e 0.0016 ± 0.0007 

obtained from other experiments. 

2. 0uPl3 DECAY IN R-PARITY NON - 
CONSERVING SUPERSYMMETRY 

R-parity of supersymmetry is defined as 

Rp = {-lfB+L+2S 

(_B, L and S are the baryon, lepton numbers and 
the spin) and is assumed to be conserved . This 
allows the R-parity conserving superpotential 
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^lH,H2, (9) 

Once R-parity is not conserved one has additional 
lepton and baryon violating terms in the super- 
potential of the form 



Q^ = 



K,kUW^Dl (10) 
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One imposes a discreet symmetry, leading to 
Kjk = or X'-jf. = 0, to avoid fast proton decay. 
In our case we will allow for lepton violation 
imposing X'^'jf. — 0. The bilinear and trilinear 
R-parity violating couplings give rise to massive 
Majorana neutrinos, which in turn contribute to 
Oiyfif] decay. Here we will consider additional con- 
tributions. The effective Lagrangian describing 
this process is 0, 0, |H1 



C 



ALe=2 



2to 



3(1 + j5)e''[r]ps JpsJps 



The color singlet hadronic currents are 

J PS = S"(l -f 75)d„, Jr = S"(T^''(l + 75)da 
•^AV = ""7^(l-75)rfa 

a is a color index and cr^'' (i/2)[7^, 7"] 
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"2 = 32/(4?^) , a. = 33/(47r) 

are SU (2) l and S'[/(3)c gauge coupling constants. 

mui^ , , mg and 

are masses of the u-squark, d-squark, gluino g and 
of the lightest neutralino Xi respectively. 
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Here and B are neutral SU{2)l and 

C/(l) gauginos while H2, are higgsinos 

which are the superpartners of the two neu- 
tral Higgs boson fields and H^. The 
neutralino couplings are defined as e^^ = 
~T:i{^)Z^ + taiiOw iT^i^l^) - Q{ip)) Zg, cr^ = 
Q{ili)ta,n9wZ g. Here Q and T3 are the elec- 
tric charge and the weak isospin of the fields 
ip = u,d, e. 

The lifetime for neutrinoless double beta decay 
now becomes: 

^ = f^n, \nm A/^^ 



GoihrA^r + (VPS - VT)Mf 



'zi'^T + ^Vps)M'^'^ + vnM 
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For comparison we write again the amphtude 
r^jvA^ivp associated with the heavy neutrino. 
The other terms come from SUSY, the first two 
correspond to the two nucleon mode, while the 
third term Ai"^^ corresponds to the pion mode. 
The pion mode is dominant PP, |H1, recently con- 
firmed with a different method 0. One gets: 
a^'^ = -0.044 , a^'' = 0.2 (Elementary particle 
treatment jSI) 

= -0.012 , a^'' = 0.15 (Quark model 
Using these ME one obtains from the data the 
limits: 

1 1 
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By employing standard techniques and working 
in the allowed SUSY parameter space one can ob- 
tain the masses and couplings involving the inter- 
mediate SUSY particles (gluinos, neutralinos, s- 
quarks and s-leptons) and their couplings, which 
enter the expression for rj^^^y. Thus from the 
limits on this lepton violating parameter, one can 
constrain and obtain limits on the unknown R- 
parity violating parameter Xm- 
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Figure 9. The R-parity violating parameter Xm 
as a function of the SUSY scale associated with 
the nucleon and the pion modes (Taken from 
Gozdz et al ^71). 
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3. BRANE WORLD 

It must have become clear from the above dis- 
cussion that Qv[3(3 decay has something to say 
in any model, which predicts massive Majorana 



Figure 10. The same as in Fig. |3|in the case of 
the pion mode for gluino {g) and neutralino (x). 
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0.05 


0.009 


0.016 


0.012 


0.008 


0.004 


0.1 


0.052 


0.061 
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0.025 


0.2 


0.096 


0.109 


0.114 


0.094 


0.058 


0.3 


0.123 


0.136 


0.143 


0.121 


0.082 


1 


0.271 


0.280 


0.287 


0.269 


0.241 


10 


0.493 
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Table 2 

Numerical estimates of (m) for different nuclei in 
a 5- dimensional brane-shifted model. (For choice 
of parameters and further details see Ji21f 



neutrinos. We cannot possibly discuss all such 
models. In models involving extra dimensions, 
however, one has the special feature that they 
predict Dirac neutrinos. Furhermore such mod- 
els predict sterile neutrinos with a mass in the 
keV-MeV region. In such cases the particle and 
nuclear physics get entangled. So the effective 
neutrino mass extracted may depend on the nu- 
cleus. 

4. NUCLEAR MATRIX ELEMENTS 

The discussion of the nuclear matrix elements 
is outside the scope of this talk. So we will limit 
ourselves to some brief comments. The first step 
towards their evaluation consists in deriving the 
effective transition operator. Particle physics dic- 
tates the structure of the operator at the quark 
level. The next step is going from the quark to 
the nucleon level. This is quite straightforward 
in the case of light neutrinos. Since, however, 
the momenta of the intermediate neutrinos are, 
by nuclear standards, quite high, w 70MeV/c, 
one must consider corrections to the structure 
of the nucleon current ^ (induced pseudoscalar 
etc), which tend to decrease the nuclear 
by 20 — 30%. This causes a corresponding in- 



crease in the extracted neutrino mass, regardless 
of the nuclear model. For heavy intermediate par- 
ticles the situation is worse. The effective oper- 
ator is short ranged. So, if one considers only 
nucleons, the nucleon must not be point like, but 
as a bag of quarks or with a suitable form fac- 
tor. Other mechanisms must also be considered, 
e.g. 01/(3(3 decay induced by pions in flight be- 
tween nucleons (pion mode). This mechanism 
tends to dominate in the case of SUSY contri- 
butions. The next step is to obtain the needed 
nuclear wave functions. This is a formidable task, 
since the nuclear systems undergoing Qv(3(3 decay 
have complicated structure. Furthermore the ob- 
tained nuclear matrix elements are quite small 
compared to a canonical value (sum rule) like 
\ME{i — > /)p (the summation is over all fi- 
nal states, not just those energetically allowed). 
So normally small effects maybe important here. 
The basic methods employed in the evaluation of 
these matrix elements Q], |5j are: 

1. The Quasiparticle Random Phase Approxi- 
mation (QRPA). One builds the needed in- 
termediate states as RPA excitations on a 
quasiparticle vacuum (one for each ground 
state). One can include a large number 
of single particle states, since the number 
of the resulting configurations is manage- 
able. This method has been applied in al- 
most all systems. Various refinements have 
been incorporated (proton neutron-pairing, 
Pauli principle corrections, gpp renormaliza- 
tion factors etc). The thus obtained matrix 
elements (ME) by various authors still have 
a spread larger than a factor of 2. 

2. The Large Basis Shell Model. With mod- 
ern computers it has been possible to use 
the nuclear shell model in obtaining the nu- 
clear wave function. One uses a single par- 
ticle model space more restricted than that 
employed in the QRPA, since soon the num- 
ber of resulting configurations becomes pro- 
hibitively large. With increasing interest in 
Ovf3(3 decay we expect the nuclear theory 
groups will improve on this in the remain- 
der of the decade, forming consortia, if nec- 
essary. 
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We thus hope that by the time the ongoing and 
planned experiments give the anticipated resuhs 
the nuclear matrix elements will be known with 
an accuracy better than a factor of 2. The re- 
liability of these matrix elements can be tested, 
once the experimental results become available, 
by taking ratios, which, assuming that one mech- 
anism dominates, are independent of the lepton 
violating parameters (ratios of lifetimes of differ- 
ent nuclei or ratios of transitions to different final 
states of the same nucleus). 

5. RESULTS 

decay is our best hope to set the scale for 
neutrino mass. Based on the available experimen- 
tal results and using the range of uncertainty in 
the nuclear matrix elements the following limits 
emerge: 

<m^> <0.45-1.4ey \'^Ge],Q.7-l.beV[^'^^Te] 

A < (0.8-2.1)xl0-6 [^^Ge],(2.5--7.3)xl0-^ P^^e,] 

7? < (0.4-1.8)xlO-*[^SG'e], (0.9-3. 7)xl0-* [^^^Te] 

Other mechanisms may dominate Qvj3[3 decay. 
One popular scenario is to be mediated not by 
neutrinos, but by SUSY particles. From rjsusY 
for each target in the order of table |21 we extract 
the value of A^: 

(1.0, 1.8, 1.8, 0.8, 1.8, 0.9, 2.3) x 10"^ 

We should not forget, however, that these values 
depend on the SUSY parameter space, in partic- 
ular the SUSY scale [201, see Figs El and CHI 

Even if the neutrinos do not dominate in Qvl3(3 
decay, its observation will demonstrate that they 
are massive Majorana particles. 
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